Non-exponential one-body loss in a Bose-Einstein condensate 



O 



S.Knoop, J.S.Borbely, R. van Rooij, and W. Vasscn 
LaserLaB Vrije Umversiteit, De Boelelaan 1081, 1081 HV Amsterdam, the Netherlands 

(Dated: February 1, 2012) 

We have studied the decay of a Bose-Einstein condensate of metastable helium atoms in an optical 
dipole trap. In the regime where two- and three-body losses can be neglected we show that the Bose- 
Einstein condensate and the thermal cloud show fundamentally different decay characteristics. The 
total number of atoms decays exponentially with time constant r; however, the thermal cloud decays 
exponentially with time constant |r and the condensate decays much faster, and non-exponentially. 
We show that this behaviour, which should be present for all BECs in thermal equilibrium with 
a considerable thermal fraction, is due to a transfer of atoms from the condensate to the thermal 
cloud during its decay. 

PACS numbers: 03.75.Hh, 34.50.Cx, 67.85.-d 
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I. INTRODUCTION 

Atomic gases can be cooled and trapped to ultracold 
temperatures and densities where Bose-Einstein conden- 
sation (EEC) and Fermi degeneracy can be reached. 
These trapped gases decay in several ways due to 
one-body, two-body and three-body collisions. Two- 
and three-body losses are density dependent and have 
been extensively studied [l|-Q and applied in work on 
atom-atom correlations and on universal few-body 
physics 0-[il- 

One-body loss is generally considered trivial as it usu- 
ally results from collisions between trapped atoms and 
room-temperature atoms from the background gas in 
the ultrahigh vacuum chamber that contains the trapped 
atoms. A background pressure of '^10^^^ mbar typically 
leads to a lifetime of atoms in the trap of ~^100 s. Other 
causes of one-body loss are scattering by off-resonant 
light in a dipole trap and on-resonance excitation by 
stray laser cooling light. All these effects cause expo- 
nential decay of the trapped cloud with a time constant 
T, typically in the range of 1 to 100 s. Experiments in 
ultracold atomic physics measure this time constant by 
monitoring the number of trapped atoms as a function of 
time, a procedure commonly performed using absorption 
imaging of the cloud on a CCD camera. Due to the high 
densities in a condensate, two- and three-body losses are 
most important in the first instances of the decay. In 
studies of decay, it is generally assumed that the decay 
becomes exponential for long enough times, where two- 
and three-body losses have become negligible due to the 
low density that is then reached. 

Zin et at. however, showed theoretically that the 
decay of a BEC is expected to be non-exponential when 
it is in thermal equilibrium with a substantial thermal 
cloud. The origin of this effect stems from the transfer 
of atoms from the condensate to the thermal cloud dur- 
ing the decay. This atomic transfer occurs only when 
thcrmalization is fast compared to the change of ther- 
modynamic variables during the decay of the trapped 
cloud. When two-body and three-body decay can be 
neglected the model predicts that for an overall expo- 



nential decay time r, the condensate decays faster and 
non-exponentially, while the thermal cloud decays expo- 
nentially with a larger time constant, |r. To our knowl- 
edge, this has not been demonstrated experimentally. 



Enhancement in the decay of a BEC in the presence 
of a thermal cloud was observed by Tychkov et at. 
In that experiment the decay of a large (> 10® atoms) 
condensate of helium atoms in the metastable 2 "^Si state 
(■*He*, lifetime 8000 s) was monitored with and without a 
thermal cloud containing approximately the same num- 
ber of atoms. That study, performed in a magnetic trap 
with atoms in the m=+l magnetic substate, revealed 
that the condensate in the presence of a thermal cloud 
decayed much faster with the cloud than without. How- 
ever, in that study the decay was studied in the presence 
of large two- and three-body losses possibly obscuring an 
effect of atomic transfer; the enhanced decay of the con- 
densate could also be understood from two- and three- 
body inelastic collisions between condensate and thermal 
atoms. Also, two- and three-body loss rate constants 
were not known accurately and both processes were ex- 
pected to contribute about equally to the decay [ll| . 
Therefore the rate constants had to be determined from 
a fit. 



Recently, we have accurately determined these two- 
and three-body loss rate constants over a range of mag- 
netic field values [l^l ■ We transferred ^He* atoms from a 
magnetic trap into an optical dipole trap and measured, 
both for atoms in the m=+l and m=-l state, the two- 
and three-body loss rate constants as a function of an 
applied magnetic field and with a (quasi-)pure BEC. In 
this paper we extend the previous experiment to par- 
tially condensed clouds with a large thermal fraction in 
order to investigate atomic transfer. To reduce two- and 
three-body losses we study long time scales and work at 
small magnetic field using m=-l atoms that, in contrast 
to m=-l-l atoms, only show three-body losses. 
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II. THEORY 



A. Theoy of trap loss 



The time evolution of the density of a trapped atomic 
gas can be described as 



n = —n/r — K2L2n'^ — n-^L^v? 



(1) 



The first term on the right takes into account one-body 
loss, which causes exponential decay with a time constant 
T. L2 and L3 are the rate coefficients for two- and three- 
body loss, respectively, and are defined such that they 
explicitly include the loss of two and three atoms per loss 
event. The constants in front of L2 and L3 are K2 = 1/2! 
and = 1 /3' for ^ BEC, while K2 = '^3 = 1 for a thermal 
gas[il. 

The rate coefficients are obtained by measuring the 
number of trapped atoms after a variable hold time. The 
analysis of this data requires integration of Eq. [1] over 
space. Because a BEC and a thermal cloud have differ- 
ent density distributions, extracting two- and three-body 
loss rate coefficients from partially condensed samples be- 
comes very difficult, and experiments usually focus on 
either a pure BEC or a thermal sample. 



B. Atomic transfer model 

In Sec. 2 of their paper, Zin et al. [l0| have derived a 
simple model for the decay of the thermal and BEC com- 
ponents of a partially condensed cloud in thermal equilib- 
rium below the temperature threshold for Bose-Einstein 
condensation. The essential assumption of the model is 
that thermal equilibrium holds during one-body decay of 
a condensate. Although two- and three-body losses may 
be incorporated in the model, this complicates the discus- 
sion and therefore we here ensure those to be negligible 
compared to the one-body losses by using low enough 
densities. If an atom is removed from the thermal cloud 
by, for instance, a collision with a background atom, then 
there is a place free in the otherwise saturated thermal 
distribution, which can be filled by a BEC atom, thus 
maintaining the size of the thermal cloud at the expense 
of the BEC. This transfer of atoms from the condensate 
to the thermal cloud enhances the decay of the conden- 
sate and increases the lifetime of the thermal cloud. 

For the whole cloud, containing N atoms, decay is ex- 
ponential with time constant r. The model then predicts 
that the thermal cloud will also decay exponentially, how- 
ever with a larger time constant, while the condensate is 
expected to decay non-exponentially [lOj : 



N = 
Nt = 

Nc = 



-N 



Nt 

At 

-\{nc + \nt 



(2) 
(3) 
(4) 



Here Nc and Nt are the number of condensed atoms 
and number of thermal atoms, respectively, and N = 
Nc + Nt- The equations show that the thermal cloud 
is expected to decay exponentially with a time constant 
r' = |t, independent of Nc^ which is assumed to be 
nonzero. For the BEC a non-exponential decay is ex- 
pected if an appreciable amount of thermal atoms is 
present. The model is valid up to the point that there 
are no BEC atoms left or thermal equilibrium cannot be 
assumed anymore. Starting with a partially condensed 
cloud finally leads to the complete depletion of the BEC; 
for the pure thermal cloud, decay then proceeds with time 
constant r' = r. Eqs. [3]and|4]are valid when the energy 
of a condensate atom (which is, apart from a small mean- 
field contribution, equal to the ground state energy eg of 
the harmonic trap potential) can be neglected compared 
to the average energy of a thermal atom {Et/Nt) [13 ■ 



cqNt/Et < 1. 



(5) 



In order to observe atomic transfer experimentally, in- 
elastic collisions within the cloud should not play a 
role, but still a high enough elastic collision rate is nec- 
essary to reach fast thermalization. The thermaliza- 
tion rate is given by 7th ~ 7coii/2.7 [l^l, with collision 
rate 7 coii = ^^av^c rd. Here, is the average density, 
V = ^JlQksT /Tim is the mean relative thermal velocity, 
and (Tci = Stto^ the elastic cross-section at ultralow tem- 



peratures, where a=142.0(0.1)ao for He* |15[. The large 
scattering length ensures fast thermalization down to a 
density of ria 



10^^ cm which is reached after 40 s 



in our dipole trap (see Sec. IIIip . 



III. EXPERIMENT 

We have studied one-body loss in a BEC of "'He* atoms 
with a considerable thermal fraction. The experimental 
setup and measurement procedure have been described 
earlier (l2l . [l6j , here we summarize only the most essential 
parts. A BEC of about 10^ atoms is prepared in a crossed 
optical dipole trap at a wavelength of 1557 nm, and all 
atoms are transferred from the m=+l magnetic substate 
to the m=-l magnetic substate by an RF sweep. We 
measure the remaining number of trapped atoms after a 
variable hold time by turning off the trap, causing the 
atoms to fall and be detected by a micro-channel plate 
(MCP) detector, which is located 17 cm below the trap 
center and gives rise to a time-of-fiight of approximately 
186 ms. From a bimodal fit to the MCP signal, the BEC 
and thermal fraction are extracted as well as the temper- 
ature and the chemical potential of the trapped gas [l6| . 

In Fig. [1] we show typical time evolutions of the par- 
tially condensed cloud, where the total, thermal and BEC 
atom number are logarithmically plotted, normalized to 
the total number of atoms at t ~ 0. Fig. [TJa) shows the 
decay for a sample with an initial BEC fraction of 80% 
(with a thermal fraction at a temperature of 0.22(3) /iK), 



3 



' ' ^ ^^^^j (a) 




10 20 30 40 

10 20 30 40 




10 20 30 40 

hold time (s) 



FIG. 1. (Color online) Lifetime measurements of an ultra- 
cold *He* gas in the m=-l state at a magnetic field of fO G, 
showing the results of bimodal fits of the MCP signal of the 
remaining ^He* atoms, where the insets show the fit temper- 
ature. The trap frequencies and the average temperature are 
(a) 2-K X 194(1) Hz and 0.22(3) ^JK and (b) 2-k x 245(1) Hz 
and 0.36(5) fxK, respectively. The solid lines and the shaded 
area are obtained from fits to the atomic transfer model, as 
described in the text. 



vifhile in Fig. [TJb) the BEC fraction is 50% (at a temper- 
ature of 0.36(5) /iK). The experiments are performed in 
a dipole trap with a geometrical mean of the trapping 
frequency of 27r x 194(1) Hz and 27r x 245(1) Hz, re- 
spectively. Eq. [5] is easily fulfilled throughout the decay 
since eoNt/Et < 0.03. We also observe that the tem- 
perature changes very little during the decay (see the 
insets in Fig. [1]). At short hold times, the loss is domi- 
nated by three-body recombination with a rate constant 
of Ls = 6.5(0.4)stat(0.6)sys x 10-^7 cm^s"! [H. Two- 
body loss by Penning ionization is strongly s upp ressed 
for the spin-stretched states m=+l and m=-l [13], while 
for m=-l the spin-dipole interaction also does not con- 
tribute to two-body losses because the atom is in the 



lowest spin state [l^. Fig. [T] clearly shows that, for hold 
times longer than 10 s, the loss of the total atom number 
becomes exponential, indicating the regime of one-body 
loss. In both measurements the initial density of the 
BEC is 3 X lO^'^ cm~^. The thermal cloud initially has 
a density riav ~ 2 x 10^^ cm~^, which after 40 s has de- 
creased to riav « 1 X 10^^ cm~'^. The thcrmalization rate 
decreases in that time from 7th « 50 s~^ to 7th ~ 3 s~^, 
still much larger than the one-body decay rate of ap- 
proximately 0.05 s^^, ensuring sufficiently rapid thermal 
equilibration. 

We observe a faster decay of the BEC fraction than the 
thermal fraction, eventually leading to a full depletion 
of the condensate, as expected from the atomic trans- 
fer model. Inspection of the decay of the thermal cloud 
for times longer than 10 s clearly shows that the ther- 
mal cloud decays with a larger time constant than the 
whole cloud. Fitting an exponential decay function to 
the data at longer hold times yields r = 18.0(8) s, r' = 
30.8(2.5) s and t'/t = 1.7(2) for the data of Fig. [^a), 
and T = 24.4(1.8) s, r' = 36.2(3.1) s and t'/t = 1.5(2) 
for Fig. [Ijb), where we give Icr-uncertainties. We at- 
tribute the relatively small lifetime to off-resonant scat- 
tering of the dipole trap light or resonant excitation by 
stray light, which could also explain the difference in the 
obtained lifetimes (our background pressure would limit 
the lifetime to ^100 s). We conclude that the thermal 
fraction decays significantly slower than the BEC frac- 
tion and that the measured lifetime ratio is in reasonable 
agreement with the theoretical prediction of t'/t = 4/3. 

To compare the measurements of the BEC fraction 
with the atomic transfer model we first fit an exponen- 
tial decay N{t) = iVoe"*/"^ to the total atom number for 
long hold times to obtain the overall lifetime r as well as 
A^o- Here Nq is the apparent atom number at t=0 in the 
absence of three-body loss. In the second step, we fit an 
exponential decay Nrit) = 7Vo(l - /c)e"3*/^^ to the ther- 
mal fraction, with only the apparent BEC fraction fc as a 
fit parameter. Finally, with the obtained parameters Nq, 
T and fc, we numerically solve Eq. |4] to obtain Nc'{t). 
The fit results are shown in Fig. [1] as solid lines. The 
Icr errors in the fit parameters are refiected in the band 
around the theoretical curve for Nc- We observe that 
the atomic transfer model describes the time evolution 
of the thermal and BEC fraction very well. 



IV. CONCLUSIONS 

Our data provides direct verification of the atomic 
transfer model. The relatively short lifetime t « 20 s, 
together with a large scattering length, provide optimal 
conditions to see this effect. Also, at the relatively small 
numbers of trapped atoms needed to demonstrate atomic 
transfer, MCP detection allows better fitting of the BEC 
fraction compared to absorption imaging. Atomic trans- 
fer is most directly visible in the one-body loss regime, 
but is important in the two- and three-body regime as 
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well, as was theoretically discussed in Ref. As a final 
remark we note that one should take care using a BEC for 
loss measurements. The thermal fraction may be small 
in the initial stage of BEC decay, but can dramatically 
increase after long hold times affecting the decay of the 
condensate, which will finally become non-exponential. 
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